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Abstract 

We consider a toy model of dark matter (DM) with a gauge singlet Dirac fermion that 
has contact interactions to quarks that differ for right-handed up and down quarks. This is 
motivated by the isospin- violating dark matter scenario that was proposed to reconcile reported 
hints of direct DM detection with bounds from non-observation of the signal in other experi- 
ments. We discuss how the effects of isospin violation in these couplings can be observed at the 
LHC. By studying events with large missing transverse momentum (Wt ), we show that the 
ratio of mono-photon and mono-jet events is sensitive to the difference in the absolute values of 
the couplings to the up and down quarks, while a dedicated study of di-jet plus Jf^ events can 
reveal their relative sign. Our methods have broad applicability to new physics that involves 
unequal couplings to up and down quarks. 



1 Introduction 



Once a hint of dark matter (DM) is found in events with large missing transverse momentum 
($r ) at the LHC, the next task will be to study how the DM particle couples to standard 
model (SM) particles. It can be the case that the DM particle couples to up and down quarks 
differently. Observing such an isospin violating feature of the new physics sector would give 
an important clue for the determination of the Lagrangian of the underlying theory. Also, it 
is itself a challenge in the physics of hadron colliders, requiring elaborate kinematical cuts to 
extract information from parton-level processes. 

In this context, isospin- violating dark matter (IVDM) [T] provides a scenario in which a DM 
particle does not couple identically to up and down quarks. This model was motivated by DM 
direct detection experiments. The DAMA [2] and CoGeNT [3] experiments observed signals 
that are potentially of dark matter origin. The signals are consistent with a DM particle of mass 
10 GeV scattering off nuclei with a spin-independent nucleon cross section of ~ 2 x 10~ 4 pb 
and ctjv ~ 7 x 10~ 5 pb, respectively. However, the XENON jl] and CDMS [5] experiments 
reported negative results. To accommodate the apparently inconsistent data, the authors of 
Ref. p] considered a dark matter particle that couples to up and down quarks differently. Then 
the cross section of dark matter scattering off protons differs from that off neutrons, and the 
apparent tension among dark matter direct detection experiments is attributed to the different 
proton-neutron ratios of the detector materials. 

We study the possibility of testing the scenario at the LHC by measuring the couplings 
to up and down quarks. We work in a general context and simply assume that a new stable 
particle that is a singlet under the SM gauge group (called "dark matter" in the following) 
couples to right-handed up and down quarks with different strengths and signs. At the LHC, 
the production of DM particles associated with a jet(s) or photon gives rise to events with a 
hard jet(s) or photon plus large missing transverse momentum )• Since the up and down 
quarks couple differently to the photon but identically to the gluons, the magnitudes of the 
two couplings can be measured by comparing the cross section of the mono-jet plus Ifii signal 
with that of the mono-photon plus $t signal. The relative sign of the dark matter couplings is 
more difficult to measure because it requires the measurement of the interference between the 
dark matter coupling to the up quark and that to the down quark. This is possible through 
the subprocess, ud — > ud + DM DM, but the interference effects need to be identified in the 
presence of the other dijet+DM DM subprocesses and the SM backgrounds. 

Our discussion of how isospin-violating couplings of a DM particle can be studied at the 
LHC involves several ingredients. We use MadGraph 5 [6] for calculating matrix elements and 
for generating DM signal events, Pythia 8 [7] for parton showering, and PGS 4 [8] for simulating 
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detector effects. The outline of the paper is as follows. In Section 2, we describe our toy model 
of IVDM. In Sections 3, 4 and 5, we investigate the DM signals from three channels, mono-jet 
plus j mono-photon plus $t , and di-jet plus • In Section 6, we compare the signals 
from the three channels and study signatures of isospin violation in the matter-DM couplings. 

2 Model 

We consider a toy model in which a DM particle couples to quarks with isospin-violating 
couplings. We introduce a Dirac fermion x that is a singlet under the SM gauge group. It 
couples to the up and down quarks through vector-like and axial vector-like contact interactions. 
The relevant part of the Lagrangian is 

£dm = iXl IJ, 9 f ,x-m x xx + -^{gqUl^L + QuUrj^ur + gDdjapd R )(xi li x) , (!) 

where A denotes the energy scale of the contact interaction and qu denotes the SU(2) doublet 
of left-handed up and down quarks. For gu ^ go, X is IVDM. 

At the LHC, we would like to extract the magnitudes and relative signs of the couplings, 
gu, go and gQ. In the following, we examine the possibility of measuring gu/gD by setting 
gQ = 0. Implications of gQ ^ will be discussed in the last section. 

3 Mono-jet + fi T 

In this section, we study events with a high pt mono-jet and large $t j which is the discovery 
channel for a DM particle at hadron colliders. We first focus on the case of pp collisions at 
\/s = 7 TeV and reproduce the results of Ref. [9]. We then consider s/s = 14 TeV, the target 
energy of the LHC. By optimizing the selection cuts, we estimate the DM signal for various 
values of gu/ 9d for a fixed value of g v + g 2 D and gQ = 0. 

3.1 Study for = 7 TeV 

Following Ref. [9], we adopt the veryHighPT cut of the ATLAS search [10], which yields the 
most stringent bound on the DM-quark interactions. This cut is described by 

• Require $ T > 300 GeV. 

• Highest p T jet satisfies |r/i| < 2.0 and pti > 350 GeV. 
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• Event is vetoed if the second hardest jet satisfies j^l < 4.5, and px2 > 60 GeV or 
A<j){p T2 , # T ) < 0.5. 

• Event is vetoed if there is another jet that satisfies \rj\ < 4.5 and pt > 30 GeV. 

• Event is vetoed if there is an electron that satisfies \r] e \ < 2.47 and px e > 20 GeV. 

• Event is vetoed if there is a muon that satisfies 1 77^ | < 2.4 and pt^ > 10 GeV. 

The dominant sources of background for mono-jet events that satisfy the above criteria are 
Z + jet events with Z — > vv, and W + jet events with W — > rv or W — > I v (I = e, fi) where 
the r decay products or leptons do not satisfy any of the veto criteria. 

We estimate the Z + jet background by first generating Z{— > vv) + 1, 2 jets events at 
the matrix element level [6] with precuts $t > 200 GeV and k? > 140 GeV for the jets. We 
next simulate parton showering [7] by using the A>r-jet matching scheme [11] to match the p? 
distribution of the second hardest jet simulated via the parton showering of Z + 1 jet events 
and the distribution calculated from the matrix elements of Z + 2 jets events. The matching 
scale is set at 200 GeV. We have confirmed that the cross section after the final cut does not 
change significantly when the matching scale is varied by ±60 GeV around 200 GeV. Finally, 
we perform a detector simulation [8] and the cross section for Z{— )■ vv) + jets events that pass 
the veryHighPT cut is found to be 

^(^)+jets (veryHighPT) ~ 106 fb . (2) 

In a similar manner, we generate W + 1, 2 jets events in which W decays into tv t and those in 
which W decays into ev e or fiv^, at the matrix element level with the precuts p l T > 200 GeV, 
where p l T denotes the transverse momentum for the three-momentum sum of the lepton mo- 
menta, and kf > 140 GeV for the jets. We then simulate parton showering with a matching 
scale of 200 GeV, and study the detector effects. Again, we have confirmed that the cross 
section after the final cut does not change much when the matching scale is varied around 
200 GeV by ±60 GeV. The cross sections for W(— > rv T ) + jets events and W(—t ev e /fiv^) + 
jets events that pass the veryHighPT cut are found to be 

(Tw{^Tu T )+)cts [veryHighPT) ~ 35.4 fb , (3) 
frwC-Ma/e/^O+joto {veryHighPT) ~ 22.2 fb . (4) 

In Table 1, we compare the numbers of background events estimated in Eqs. (T2JHj) with 
those estimated by the ATLAS collaboration [10] for an integrated luminosity of 1 fb -1 . We 
find satisfactory agreement. Following Ref. [9], we normalize our estimates for Z ± jets events 
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Background 


Our simulation 


ATLAS collaboration 


Normalization factor 


Z{yv) + jets 


106 


124 ±12 ±15 


1.17 


W(tu) + jets 


35.4 


36±7±8 


1.08 


W{e/jiv) + jets 


22.2 


26 ± 4 ± 3 



Table 1: The numbers of background events that pass the veryHighPT cut for 1 fb _1 LHC (at 
7 TeV) according to our simulation (see Eqs. IHH]), compared to the estimates by ATLAS [TO] . 
The last column gives the ratios of the ATLAS estimate to our estimate. 

and for W + jets events to those of the ATLAS collaboration. The normalization factors are 
also shown in Table 1. Note that the normalization factor for the Z{yv) + jets background 
(1.17) is consistent with Ref. [9], where it is estimated to be 1.2. 

The DM signal in mono-jet events arises mainly from the following subprocesses: 

9 u -> u x X , (5) 
9 d ->■ d x X • (6) 

We simulate the DM contribution to mono-jet events by following the method we have adopted 
for the background simulations. We generate pp — > XX + 1> 2 jets events at the matrix element 
level with the precuts Ifii > 200 GeV and > 140 GeV for the jets. Parton showering is 
simulated with a 200 GeV matching scale, and detector effects are simulated. We have also 
confirmed for the DM signal that the final results are not significantly affected by varying the 
matching scale around 200 GeV. Finally, we rescale our signal by a factor of 1.17; see Table 1. 

In order to compare our results with those of Ref. [9], we estimate the DM signal with 
A = 400 GeV, gu = 9d = 9q = ^ an d m x = 10 GeV. We show the distribution in Fig. [TJ 
All the curves are in good agreement with those in Fig. 2 of Ref. [9]. The difference in the 
300-350 GeV bin is likely due to the different modeling of parton showers. 

To compare the lower bound on A from our analysis with that from Ref. [9], we derive the 
observed bound for m x = 10 GeV and gu = Qd = Qq = 1- It should be noted that only when all 
three couplings of the effective Lagrangian in Eq. ([[]) are equal, do we have purely vector-like 
couplings for both up and down quarks, as assumed in Ref. [9]. Isospin violation with purely 
vector-like couplings is not compatible with electroweak symmetry. To place 90% C.L. lower 
limits on A we require [9] 



2 = (N obs - N DM (A) - N SM ) 2 m 
N DM (A) + N SM + a 2 SM 



where N b s ,NDM and Nsm are the number of observed events, the expected number of signal 
events, and the SM expectation, respectively, and as\i is the uncertainty in N SM . We adopt 
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Figure 1: The $t distribution of the backgrounds and the DM signal for 1 fb _1 at ATLAS. 
The curves from bottom to top successively include the W{e/jiu) + jets background, the 
W(tv) + jets background, the Z[yv) + jets background, and the DM signal with A = 400 GeV 
and m x = 10 GeV. i.e., the curve labelled "DM" corresponds to the sum of all the backgrounds 
and the DM signal. 

N bs — 167 and N$m =t &sm — 193 ± 25 as obtained by the ATLAS collaboration in 1 fb" 1 of 
data. After rescaling, we find the number of signal events for A = 400 GeV to be Ndm = 319, 
which implies: Ndm {A) = 319 x (400GeV/A) 4 . Finally, using the criterion Eq. ([7]), we obtain 
the 90% C.L. bound, 

A > 783 GeV . (8) 

On the other hand, Fig. 4 of Ref. [9] shows separate bounds on A for DM-up quark and DM- 
down quark interactions, which are both purely vector-like. These bounds can be translated 
into a bound on A for gu — go = Qq — 1: 

A > (700 4 + 575 4 ) 1/4 GeV = 769 GeV , (9) 

which is consistent with our estimate in Eq. (JSJ). It is clear that our estimate of the number of 
signal events is consistent with that of Ref. [9]. 

For gy — 1 and gjj = gQ = 0, we derive the following bound on A using the data and the 
SM background estimates of the ATLAS collaboration [10] : 

A > 599 GeV . (10) 
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3.2 Study for ^/s = 14 TeV 

In this subsection, we study the DM contribution and the SM background to the mono-jet 
cross section at the LHC with yt~s = 14 TeV. We introduce the following cut, which we refer to 
as MonoJUTeV: 

• Require $ T > 800 GeV. 

• Highest pt jet satisfies \rji\ < 2.0 and pn > 700 GeV. 

• Event is vetoed if the second hardest jet satisfies |?7 2 | < 4.5, and px2 > 120 GeV or 
A(j)(p T2 , E T ) < 0.5. 

• Event is vetoed if there is another jet that satisfies \rj\ < 4.5 and pt > 60 GeV. 

• Event is vetoed if there is an electron that satisfies \r] e \ < 2.47 and pt £ > 20 GeV. 

• Event is vetoed if there is a muon that satisfies 1 77^ | < 2.4 and Pt^ > 10 GeV. 

In comparison to the ATLAS veryHighPT selection cut, we keep the lepton veto conditions 
and scale all the jet pt cuts by a factor of 2, while the $t cut is chosen to reduce the Zivv) 
+ jet background more efficiently in order to increase the signal significance for smaller DM 
signal cross sections. 

We estimate the Z{yv) + jet background in a manner similar to that for -y/s = 7 TeV, by 
changing the precuts to $r > 600 GeV and kr > 240 GeV for the jets, and the matching scale 
to 340 GeV. After testing the stability of our estimate under variation of the precuts and the 
matching scale, we find that the Z + jets background cross section is given by 

a z+jets (MonoJUTeV) = 12.8 fb . (11) 

Similarly, we estimate the cross sections for the W{rv) + jet and W{lv) + jet backgrounds with 
the precuts p l T > 600 GeV and kr > 240 GeV for the jets, and a matching scale of 340 GeV. 
We find 

a W (Tu)+jets{MonoJUTeV) = 2.4 fb , (12) 
a w{lu)+jets (MonoJUTeV) = 2.4 fb . (13) 

To estimate the DM signal contribution, we set {gQ, gu, 9d) = (0, cos (ft, sin 0) so that the 
magnitudes of the couplings are given by 1/A 2 , and the isospin violation is parametrized by the 
phase — 7r/2 < 4> < ir/2. The up and down quark couplings have the same sign for < (f> < tt/2, 
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m x 


= 


= ±7r/4 


= ±7T/2 


10 GeV 


15.9 


11.0 


6.37 


300 GeV 


15.1 


10.5 


5.91 


500 GeV 


13.0 


8.89 


4.98 



Table 2: DM signal cross sections (in units of (800GeV/A) 4 fb) for yfs = 14 TeV with the 
MonoJIATeV cut, for various values of the DM mass m x and the ratio of the DM couplings 
parametrized by = tan -1 (go/gu)- 

and opposite signs for — n/2 < < 0. We repeat the simulation steps of the y/s = 7 TeV case, 
by changing the precuts to $t > 600 GeV and kx > 240 GeV for the jets, and the matching 
scale to 340 GeV. We calculate the cross sections for various values of the DM mass m x and 
= tan -1 ^ o/gu)] some of these are listed in Table 2 in units of (800 GeV/A) 4 fb. The mono- 
jet cross section does not distinguish between the signs of the gu and go couplings, and hence 
±0 give the same prediction. 

4 Mono-photon + fa 

The mono-jet cross section alone does not contain information on 0. Therefore additional 
signals are needed, and we study the mono-photon plus $t signal in this section. Since the 
photon coupling to the up quark is twice that to the down quark, we expect a strong dependence 
of the DM signal on |0|. We first focus on the \fs = 7 TeV case and reproduce the results of 
Ref. j9]. We then proceed to the y/s = 14 TeV case, and look for an appropriate selection cut 
at high energies and at high integrated luminosity. 

4.1 Study for = 7 TeV 

Following Refs. (9j [12], we implement the following selection cut, which we call MonoG7TeV: 

• Require a photon with |t/ 7 | < 1.44 and pt^ > 95 GeV. 

• Require $ T > 80 GeV. 

• Event is vetoed if there is a jet with \r]j\ < 3.0 and pxj > 20 GeV. 

• Event is vetoed if there is an isolated lepton with p Ti > 10 GeV and AR(l,j) = 
V(Vi-^) 2 + (<Pi-<Pj) 2 > 0.04. 
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The leading background arises from Ziyv) + 7 events. In addition, Zivv) + jet events with 
the jet misidentified as a photon, and W(ev) events with the electron misidentified as a photon 
also contribute to the background. 

We simulate the Zivv) + 7 background as follows. We first generate Zivv) + 7 and Ziyv) 
+ 7 + 1 jet events at the matrix element level with the precuts py 7 > 60 GeV for the photon 
and kr > 60 GeV for the jet, and simulate parton showering by using fcy-jet matching with 
a matching scale of 84 GeV. After simulating the detector effects with PGS and imposing the 
selection cut, we obtain 

a Z7 (MonoG7TeV) = 48.0 fb . (14) 

This gives 54.7 events in 1.14 fb -1 of data, while the CMS collaboration estimates 36.4 events [T2] . 
We therefore multiply the DM signal (which has a similar event topology as the Ziyv) + 7 
background) by the ratio 36.4/54.7 = 0.67. This rescaling reflects the additional criteria for 
photon indentification specific to the CMS detector. Note that the factor of 0.67 is in agreement 
with the analysis in Ref. [9], which finds 0.71. 

The DM signal arises mainly from the following subprocesses: 

u u -»■ 7 x X , (15) 
d d -»■ 7 x X • (16) 

Setting m x = 10 GeV and gu = 9d = 9q = 1, we generate XX + 7 an d XX + 7 + 1 jet 
events at the matrix element level with the same precuts as for the background, which are then 
processed to the parton showering simulation with the same matching scale and finally the 
detector simulation is done using PGS. In this way, we find the DM signal cross section before 
the rescaling to be 

a DMj (MonoG7TeV) = 63.3 x (400GeV/A) 4 fb . (17) 

To compare our analysis with that of Ref. [9], we estimate the 90% C.L. lower bound on A by 
using the criterion of Eq. (0). For Nsm and <tsm, we adopt the numbers estimated by the CMS 
collaboration, which are Nsm — 67.3 and <jsm — 8.4 in 1.14 fb -1 [12]. The CMS collaboration 
reports N t, s = 80 [12]. Njjm is estimated by multiplying Eq. (fTTj) by the normalization factor, 
0.67. We obtain the 90% C.L. bound, 

A > 422 GeV . (18) 

On the other hand, Fig. 8 of Ref. [9] shows separate bounds on DM-up quark and DM-down 
quark interactions. These bounds can be translated into a bound on A for the case, gu = 9d = 
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Cut 


0"Z 7 




&DM~ f / y/fZ-y + CTDA/7 


MonoGUTeV-a 


35.8 


5.49 


0.85 


MonoGUTeV-b 


11.6 


3.34 


0.87 


MonoGUTeV-c 


4.51 


2.13 


0.83 



Table 3: Cross sections (in fb) and significance factors (fb 1 / 2 ) for the DM signal with 
m x = 10 GeV, A = 800 GeV and = 7r/4. The three selection cuts for mono-photon events 
are compared. 

9q = 1: 

A > (400 4 + 240 4 ) 1/4 GeV = 412 GeV . (19) 
This is consistent with our result in Eq. ( TT8|) . 

4.2 Study for ^ = 14 TeV 

We introduce the following selection cut for the LHC with y/s = 14 TeV, which we name 
MonoGUTeV-a: 

• Require a photon with |r/ 7 | < 1.44 and pr 7 > 140 GeV. 

• Require $ T > 140 GeV. 

• Event is vetoed if there is a jet with \rjj\ < 3.0 and p?j > 40 GeV. 

• Event is vetoed if there is an isolated lepton with p^i > 10 GeV and Ai?(Z, 7) > 0.04. 

We also introduce the selection cut M onoGlAT eV -b which requires (pt 7 , $t) > (200, 200) GeV, 
and the cut MonoG14TeV-c which requires (pry, $t) > (260, 260) GeV, and for which all the 
other requirements are the same as for MonoGl4TeV-a. 

We estimate the Zivv) + 7 background in a similar manner to the case of \fs = 7 TeV, 
by changing the precuts to p Tl > 100 GeV for the photon and k T > 100 GeV for the jet, and 
the matching scale to 140 GeV. We estimate the DM signal contribution with the same precuts 
and matching scale as for the Z{yv) + 7 background estimation for = 14 TeV. 

Table 3 shows the cross sections and significance factors (S/ V S + B) for the DM signal 
with m x = 10 GeV, A = 800 GeV and <p = n/4, where the three selection cuts are imposed. 
The numbers are in units of fb for the cross sections and fb 1//2 for the significance factor. We 
find that the significance factor does not depend much on the choice of the selection cut. This 
is because the background cross section decreases by a factor of 8 and the signal also decreases 
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Figure 2: The Z + 7 background cross section and the DM signal cross section for A = 800 GeV 
that satisfy the MonoGlATeV-a criteria, as a function of |0| for the couplings {gQ, gu, 9d) = 
(0, cos0, sin 0). The horizontal solid line corresponds to the background cross section. The 
curves correspond to the background + signal cross sections for three values of m x . 

by 2.6 ~ v8 when we change the cut from MonoGlATeV-a to MonoGlATeV-c. For our study, 
we adopt the selection cut MonoG14TeV-a which yields the largest signal cross section. This 
is because, with more events, it is easier to discriminate the signal from the background from 
the difference in their and distributions as well as from the correlation between 7 and 
$t momenta, although we do not perform such an analysis in this paper. 

In Fig. [21 we show the cross sections for the Z + 7 background and the DM signal for 
A = 800 GeV and (gQ, gu, go) = (0, cos0, sin0) that satisfy the cut MonoG14TeV-a. The 
three curves correspond to m x = 10, 300, 500 GeV. The mono-photon cross section does not 
depend on the signs of the gu and go couplings, and hence ±0 give the same prediction. 

In Table 4, we present the DM signal cross sections in units of (800 GeV/A) 4 fb for three rep- 
resentative couplings, = (gu = 1), = ±7t/4 (\gu\ = \go\ = l/v^2) and — 7r /2 (gD = !)• 
Roughly speaking, the = n/2 (gu, = 1) case gives about a factor of 7 smaller DM signal cross 
section than the = (gu = 1) case, while the cross section in the = ±7r/4 (\gu\ = \go\ = 
1/V2) case is about half of that in the = case, for all m x between 10 GeV and 500 GeV. 
This reflects the combined effect of the QED coupling ratio (Qd/Qu) 2 = 1/4 and the ratio of 
the down quark to up quark parton distribution functions in pp collisions. 
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m x 


= 


(f) = ivr/4 


(f) = ±vr/2 


10 GeV 


9.69 


5.49 


1.38 


300 GeV 


8.69 


4.96 


1.22 


500 GeV 


6.98 


3.98 


0.945 




0.0 1 1 1 1 1 1 1 

15 30 45 60 75 90 

101 [degree] 



Figure 3: Ratio of the DM signal cross section in the mono-photon channel with the 
MonoG\4TeV-a cut to the cross section in the mono-jet channel with the MonoJIATeV cut 
for three values of m. 
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In Fig. |3l we show the ratio of the DM signal cross sections in the mono-photon channel with 
the M onoGlAT eV -a cut to the cross section in the mono-jet channel with the MonoJ14TeV 
cut. We find that this ratio is very sensitive to the coupling ratio \gol9v\ — tan|0|, but 
insensitive to the DM mass m x . It is independent of A and \/ gfj + g 2 D . 

Experimentally, the cross section ratio, o- 1 DM /a 3 DM) is determined by 

crJ M (monophoton cut) (^L — iV^ M )(monophoton cut) 



(20) 



cr^ M (monojet cut) {^Ibs ~~ ^sAi)( mono i e ^ cu ^) 

where N^ bs and N^ bs are the observed number of mono-photon and mono-jet events, respectively, 
and Ng M and N J SM are the SM expectations for mono-photon background events and mono-jet 
background events, respectively. Since the cross section ratio has a common value for a wide 
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range of m x , we can determine the value of |0| by comparing the observed value of the cross 
section ratio with Fig. [3j 

As an illustration, we estimate the integrated luminosity L needed to measure the cross 
section ratio with 10% accuracy close to the point \<ft\ = n/4 for A = 800 GeV and m x = 10 GeV. 
We assume that the statistical uncertainty of the number of events, N, follows AiV = \^N, and 
ignore systematic uncertainties. To estimate N^ bs and Ng M , we rescale the cross sections by the 
normalization factor 0.67 obtained in Section 4.1. The statistical uncertainty of udm^I '&DMj is 
given by 

a DM 

Inserting NJ M = L x 0.67 x 35.8 fb and N^ bs - N] M = L x 0.67 x 5.49 fb from Table 3, 
N 3 SM = L x (12.8 + 2.4 + 2.4) fb from Eqs. (JUS), and N j ohs - N j SM = L x 11.0 fb from Table 2, 
we find that in order to obtain a 10% measurement of the monophoton-to- mono jet cross section 
ratio, L ~ 420 fb -1 is needed. 




5 Di-jet + £ T 

Since the mono-jet and mono-photon cross sections depend only on the absolute values of gu 
and go, a measurement of the relative sign of gu and go is not possible from these channels. 
To determine the relative sign, we now focus on events with two hard jets and large ^ . 
This channel can be sensitive to the relative sign of the gu and gu, couplings because in the 
subprocess, 

u R d R u R d R xx , (22) 

the amplitudes where XX are emitted from the up quark interfere with the amplitudes where 
XX are emitted from the down quark, as can be seen from the Feynman diagrams in Fig. HI 
The interference term is directly proportional to gugo so that the cross section depends on the 
sign of g v g D . 

5.1 Extracting the Interference 

In this subsection we discuss kinematic conditions that enhance interference effects in the 
diagrams of Fig. HI 
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Figure 4: Feynman diagrams for the subprocess, urcIr — > UrcLrxX- The diagrams on the left 
are proportional to gu, and the ones on the right are proportional to go- 

5.1.1 Suppressing contributions from quark-gluon collisions 

First, it is necessary to suppress large contributions from quark-gluon collisions, 

q 9 -> q 9 Z(w) , (23) 
q 9 -> q 9 XX , (q = u, d). (24) 

While the subprocesses in Eq. (1231) give the dominant background for the DM plus di-jet 
production events, the subprocesses in Eq. (124"1) should also be suppressed because their cross 
sections do not depend on gugD- The contributions from quark-gluon collisions can be reduced, 
compared to the signal subprocess in Eq. (1221) from quark-quark collisions by imposing a very 
large cut on the "di-jet cluster transverse mass" [T3] : 

M T (jj;# T ) = yf + \p 1T + fcr\ 2 + #t , (25) 

where M\i is the invariant mass of the two hardest jets, and pit +P2T is the vector sum of the 
transverse momenta of the two hardest jets. Since the parton center-of-mass energy (\/S) is 
always larger than M T (jj; $r)) by imposing a cut on M T (jj; we can select those events 
that satisfy 

VI > M T (jj;# T ) > (M r (jj;^ T )r f . (26) 
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By choosing a large value for (Mr(jj; Ifii)) cut , we can suppress the contributions from quark- 
gluon collisions relative to those from quark-quark collisions, because the former subprocesses 
have relatively less parton center-of-mass energy. We note that the contributions from the same 
quark collisions: 

q q q q Z{vv) , (27) 

q q -> q q xx , (q = u, d), (28) 

cannot be suppressed by a cut on Mt(jj;-^t)- We therefore treat them as irreducible back- 
grounds for a measurement of the sign of gugr>- 

5.1.2 Utilizing the generalized null radiation zone theorem 

Generally speaking, when a subprocess contains two Feynman diagrams that interefere, the 
ratio of the interference term to the cross section is maximized when the amplitudes of the 
two diagrams take the same absolute value. The null radiation zone theorem [H] provides us 
with powerful criteria for identifying such kinematic regions. The theorem states that for any 
tree- level Feynman diagram, if the ratio Qk/{Pk • q) (where Qk is the charge of an external 
particle, pk is its four-momentum, and q is the four-momentum of the emitted photon) is 
the same for all external particles and if the charge Qk's are conserved, then the sum of the 
amplitudes of all the tree-level diagrams made by adding one photon emission vertex to the 
original diagram vanishes for all helicities. Since the amplitude for each diagram does not 
vanish, the theorem implies that the amplitudes of the contributing diagrams cancel exactly, 
i.e., they interfere maximally. The original theorem applies to massless vector boson emission 
via a vector coupling to fermions. We make use of the theorem for our dark matter pair which 
couples to quarks via vector and axial vector couplings. Since the dark matter pair is not 
massless, we first generalize the theorem to emission of a massive vector current. 

We now show that the null radiation zone theorem can be generalized to emission of a 
massive vector boson, or a vector current whose invariant mass squared is time-like, q 2 > 0. 
Consider the following process in which a neutral vector current, V, is emitted: 

a + b ->■ 1 + 2 + ... +n + V . (29) 

We label the particle four-momenta and charges in the initial state by Pi, Qi (i = a,b), and 
those in the final state by pf, Qf (/ = 1, 2, n). We denote the four-momentum of the 
vector current by q. The tree-level scattering amplitude should vanish for all helicities when 
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the following conditions are satisfied: 

Qi Qf 



2pi - q — q 2 tyf ■ q + q 2 



(a common value) for all i and / , (30) 



j2q, = ^TQ. f . (3i) 

i f 

Let us focus on the four diagrams of Fig. HI We denote the four-momenta of the incoming 

quarks by k\ and k 2 , those of the outgoing quarks by p\ and p 2 and the four-momentum sum 

of the DM momenta by q. We notice that the values of 2pi ■ q + q 2 and 2p 2 ■ q + q 2 are always 

positive. The values of 2k\ ■ q — q 2 and 2k 2 ■ q — q 2 can take both signs, but their sum is always 

positive because k\ + k 2 = p\ + p 2 + q- Therefore the null radiation zone can be realized only 

when gu and go take the same sign. We thus expect that the interference in pp collisions is 

destructive if gu and go take the same sign, and is constructive otherwise. 

The theorem suggests that if the condition, 

q2 ^2 q2 ^2 

h-q~— = h-q~ — =Pi-q+ — =P2-q+~, (32) 

is satisfied, the amplitudes of the four diagrams of Fig. H] cancel completely for gu = go, while 
the sum of the amplitudes is maximally enhanced for gu = —go- Equation (132]) reduces to the 
following kinematic conditions for massless partons satisfying p\ = p\ = k\ = k\ = 0: 

\Pi\ = \P2 1 , (33) 
q 2 = 0, (34) 

where and \p 2 \ are respectively the magnitudes of the three- momenta of p\ and p 2 in the 
colliding parton center- of -mass frame. Since our DM particle is not massless, the condition 
Eq. ( 1321) cannot be satisfied and the null radiation zone does not exist in the physical region. 
However, we expect that strong destructive interference occurs for gu = go when the condition 
Eq. (152"]) is approximately satisfied. 

At hadron colliders, we cannot measure the missing mass \/q^, and cannot determine the 
colliding parton center-of-mass frame. Therefore, to enhance kinematic regions around the null 
radiation zone, we must make use of the constraints among jet transverse momenta coming 
from Eq. (133]) . In the following we examine the selection cut, 

\pri ~Pt2\ < C pti , (35) 

where p?i and pT2 are the transverse momenta of the hardest and second hardest jets, respec- 
tively, and C is a number less than unity. Although no selection cut can be applied to enhance 
kinematic regions that nearly satisfy Eq. (132]) . such regions are automatically favored in pp 
collisions with a fixed cut because they are the regions where the parton center-of-mass energy 
is minimized. 



16 



5.2 Analysis 

The "signal" in this analysis is not DM production itself, but the difference in the DM signal 
cross sections for gu /9d > and gu/ go < 0. In this analysis, therefore, we consider the value of 
o'DMjji'fi — ~' K /4)—a£)Mjj{. ( t> — as t ne signal and treat the value of crDMjj(<i> — ~ ^ /4)+ctsm 
as the background. 

We first examine selection cuts dubbed DiJ-a, DiJ-b, DiJ-c and DiJ-d. The DiJ-a cut is: 

• Require two jets with \r)\ < 4.5 and p? > 200 GeV. 

• Require # T > 300 GeV. 

• Require M T (jj;# T ) = ^ Mf 2 + \p 1T + p 2T \ 2 + tf T > 2 TeV. 

• Event is vetoed if there is another jet that satisfies \rj\ < 4.5 and pt > 100 GeV. 

• Event is vetoed if there is a jet whose three-momentum p satisfies A<f)(p, ^t) < 0.2. 

The first, second and fourth conditions define the two jets plus large $r events that we study. It 
is the third condition that enhances events from quark-quark collisions over those from quark- 
gluon collisions. The last condition is necessary to reduce those Z[yv) + jets events in which a 
quark emits the Z boson almost collinearly; the amplitude receives collinear enhancement even 
for the Z boson because Ej et /Mz 3> 1. For comparison, we consider a selection cut DiJ-b that 
requires Mt{jj] $t) > 2.5 TeV, and a selection cut DiJ-c that requires Mt(H] $t) > 3 TeV, 
with the other conditions the same as for DiJ-a. Based on the null radiation zone theorem, we 



also consider a selection cut DiJ-d that requires 

\p T1 -p T2 \ < 0.5 p TU (36) 

in addition to the conditions of DiJ-b. (pti and pT2 are the transverse momenta of the hardest 
and second hardest jets, respectively.) The cuts are summarized as follows: 

DiJ-a : M t (jj;^ t ) > 2 TeV . (37) 

DiJ-b : M T {jj;tf T ) > 2.5 TeV . (38) 

DiJ-c : M T (]j;tf T ) > 3 TeV . (39) 

DiJ-d : M T (jj;# T ) > 2.5 TeV and \p T1 - p T2 \ < 0.5 p T1 . (40) 



Our simulation of the SM background only accounts for the dominant background from 
Z{yv) + jets events. To estimate the background cross section, we generate Z{yv) + 2, 3 jets 
events at the matrix element level with the following precuts: > 200 GeV, > 160 GeV, 
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ST rr ■ ■ 


ODM]]\<P — 4 J 


/T 7-4 71 /r ■ ■ ( A} — 1 

C>DMjj\9 — 4 J 




o-DA/jj(-7r/4)-cr DA /jj(7r/4) 






DiJ-a 


49.8 


9.246 


9.595 


0.0364 


0.0453 


DiJ-b 


20.1 


4.094 


4.355 


0.0599 


0.0528 


DiJ-c 


8.34 


1.868 


1.999 


0.0655 


0.0407 


DiJ-b 


13.7 


1.847 


2.012 


0.0820 


0.0416 



Table 5: Cross sections (in fb), the dimensionless ratio of the interference effect to the DM pro- 
duction cross section, and the significance factor (fb 1//2 ) for DM production with m x = 10 GeV, 
A = 800 GeV and = ±vr/4, for four different cuts. 

AR(pj,p z ) > 0.2, V§ > 2 (2.5) TeV when using the DiJ-a cut (DiJ-b, DiJ-c and DiJ-d 
cuts), where pj is the three-momentum of any jet and pz that of the Z boson. We match the 
matrix element events with the parton shower with a matching scale of 220 GeV, and perform a 
detector simulation. We find that the cross section after the final cut does not change drastically 
when the matching scale is varied by ±40 GeV. 

The DM production cross section is estimated in a similar manner; we generate XX + 2, 3 
jets events with the same precuts as the background simulation (pz is replaced with the three- 
momentum sum of the DM momenta), and process them to parton showering and detector 
simulation with the same matching scale. 

Table 5 shows the cross sections and the significance factors (S/\/S + B) for DM production 
with m x = 10 GeV, A = 800 GeV and = ±7r/4, for the four selection cuts. Also shown is the 
ratio of the interference effect (<JDMjj( ( t ) = — 7r /4) — &DMjj{4> — + 7r /4)) to the DM production 
cross section for = — 7r/4 (o'DMjj{<j ) = ~ 7r / / ^))- Notice that in this analysis the "signal" 
cross section corresponds to the difference between the DM cross sections for <ft = ir/4 and 
<p = — tv /A. The numbers are in units of fb for the cross sections and fb 1 / 2 for the significance 
factors. From Table 5 we see that the DiJ-b cut gives the largest significance factor. We note 
in passing that the requirement Eq. (1361) of the DiJ-d cut, which is based on the null radiation 
zone theorem, does enhance the ratio of the interference effect to the DM production cross 
section, {(TDMjj(4> = ~ n /^) ~ ® DMjj (0 — + 7r /4)}/cr£>A'/jj(0 = — 7r/4). Unfortunately, the Z + 
jets background is not significantly diminished by this requirement, resulting in a significance 
factor that is smaller for the DiJ-d cut than for the DiJ-b cut. 

In Fig. [5] and Table 6, we show cross sections for the Z + jets background and the DM 
production process for A = 800 GeV and m x = 10, 300 GeV that satisfy the DiJ-b selection 
cut. The numbers in the table are in units of (800 GeV/A) 4 fb. Finally in Fig. EJ we show the 
ratio of the DM signal cross section in the di-jet channel with the DiJ-b cut to the cross section 
in the mono-jet channel with the MonoJ14TeV cut, as a function of 0. 
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Figure 5: The Z + jets background cross section and the DM production cross section 
for A = 800 GeV that satisfy the DiJ-b selection cut, as a function of for the couplings 
{9Qi 9Ui 9d) = (0, cos0, sin0). The horizontal solid line corresponds to the background cross 
section. The solid and dotted curves correspond to the background + signal cross sections for 
m x = 10 GeV and 300 GeV, respectively. For each pair of curves, the upper one corresponds 
to < and the lower one to > 0. 





= +vr/4 


= -7T/4 


10 GeV 


4.094 


4.355 


300 GeV 


3.712 


3.895 



Table 6: DM production cross sections (in units of (800GeV/A) 4 fb) for a/s = 14 TeV that 
satisfy the DiJ-b selection cut. 
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Figure 6: Ratio of the DM production cross section in the di-jet channel with the DiJ-b cut to 
the cross section in the mono-jet channel with the MonoJ14TeV cut for two values of m x . 



5.3 Discussion of the Di-jet Channel 

From Fig. [5] and Table 6, we confirm that the Feynman diagrams of the subprocess, urcIr — > 
XX u RdR, interfere destructively if gu and go have the same sign, and constructively if they 
have opposite signs. This is in accordance with our expectation in Section 5.1 based on the 
null radiation zone theorem [13]. Table 6 also indicates that the effect of the interference is 
diminished for larger DM masses. This is because for heavier DM, the physical region deviates 
more from the null radiation zone given by Eq. (132|) . 

The deviation of Fig. [6] from a left-right symmetric form is a consequence of the interference. 
Apart from the left-right asymmetry, Fig. |6] shows that the cross section ratio is smaller for 
(f) ~ ±7r/2 than for <fi ~ 0. This is due to the difference in the parton distribution functions 
of the up and down quarks; since the momentum distribution of the down quark in a proton 
leans towards a smaller momentum region compared to that of the up quark, events involving 
only down quarks have relatively less parton center-of-mass energy and are more likely to be 
rejected by the high cut on Mr(jj; -^t)- 

Once it has been experimentally established that \gu\ — \go\ (1 01 — 7r /4), and that 



m x ~ 10 GeV, the sign of gugD can be determined by measuring the following DM cross section 
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ratio and comparing it with the prediction of Fig. [6j 



n 

DM 



(Nii s -Ni? M )(DU-b) 



^dm (NL - N 3 SM ){ManoJUTeV) 



(41) 



where N^ s is the observed number of di-jet events and Njg M is the SM expectation for di-jet 
events. 

We estimate the integrated luminosity needed to determine the sign of gugD for A = 800 GeV, 
m x = 10 GeV and \<fr\ = tt/A. To distinguish <p = — 7r/4 from <fi = +7r/4 at the 2a level, the 
dijet-to- mono jet cross section ratio must be measured with 3% accuracy since Table 5 shows 
that the difference between the 2 cases is about 6%. The statistical uncertainty of (JDMjjl &DMj 
is given by 

± (°DMn\ ( °DMn \ I K bs + N> SM Kl + Ng~ 

V °dmi > V mm ) V (K bs - ^ SM ? (Kl - Ni j M y ' 

Inserting N j s j M = L x 20.1 fb and N 3 J S - N j s j M = L x 4.355 fb from Table 5, and N j SM = 
L x (12.8 + 2.4 + 2.4) fb from Eqs. (jIUH3j) and N J obs - N j SM = L x 11.0 fb from Table 2, we find 
that to achieve a 3% measurement of the dijet-to- mo no jet cross section ratio, L ~ 3000 fb -1 
is needed. 



6 Conclusions 

We explored ways of observing the isospin violating nature of a dark matter particle at the 
LHC We adopted a Dirac fermion DM particle that couples to up and down quarks through 
contact interactions with different strengths. The cross section ratio of the DM production 
associated with a hard jet and that with a photon reflects the absolute value of the ratio of the 
DM couplings to up and down quarks because they have different electric charges. We showed 
that this ratio can be measured at the LHC by observing and comparing mono-photon + ^ 
events with mono- jet + ^ events. 

The relative sign of the DM couplings to up and down quarks can be studied only by com- 
paring the cross section of DM production associated with two hard jets to that with one hard 
jet. This is because the subprocess, urcI r — > undnxXi that contributes to di-jet + ^ T events, 
contains an interference term that is proportional to gugD- in order to extract the interference 
term, we implemented a very hard cut on the di-jet cluster transverse mass, Mr(jj; -^t), which 
selects quark-quark collision events over quark-gluon collision backgrounds. We showed that 
the effect of the interference does in fact appear in the cross section of di-jet + Ifii events with 
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the sign expected from the null radiation zone theorem [13], and that it is possible to determine 
the relative sign of the couplings by observing and comparing di-jet + fa events with mono-jet 
+ fa events. 

Throughout, we focused on the case with gQ = 0. If gQ ^ 0, the DM cross section ratios in 
Figs. [3] and [61 will be modified. For Fig. [3], the angle will be replaced by 

<f>eff = tan^ 1 ( yJg 2 D + g 2 Q / yj 9u + 9q ) ■ 

On the other hand, Fig. [6] will be modified in a complicated way. Note that the ratio of the 
di-jet and mono-jet signal cross sections can be written in the following way when gQ = 0: 

Cgfj + Dgygp + Eg 2 D 

^ = + ' (43) 

where A, B, C, D, E are numerical factors with A, B, C, E > and D < as we saw in Section 
5. For gQ ^ 0, the expression above becomes: 

C{gv + g 2 Q ) + D{gug D + g Q gu + g Q g D ) + E(g 2 D + g Q ) 
A{glr + g 2 Q ) + B{gl + gl) 

Our analysis can be generalized to the entire parameter space of (gQ, gu, go)- 

Although this work is motivated by isospin violating dark matter, our methods are applica- 
ble to other new physics models for which one studies the ratio of the new physics couplings to 
up and down quarks. The di-jet channel is challenging because it is necessary to suppress large 
contributions from gluon-quark interactions. However, it is the only channel that is sensitive 
to the relative sign of the up and down quark couplings, and hence should be studied seriously 
once new physics that couples to quarks is discovered. We believe that our exploratory studies 
will turn out to be useful to probe new physics properties at the LHC. 
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